Quantum control of electron— phonon scatterings in artificial atoms 
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The phonon-induced dephasing dynamics in optically excited semiconductor quantum dots is 
studied within the frameworks of the independent Boson model and optimal control. We show that 
appropriate tailoring of laser pulses allows a complete control of the optical excitation despite the 
phonon dephasing, a finding in marked contrast to other environment couplings. 
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Ultrafast spectroscopy allows to investigate the fun- 
damental interaction mechanisms in semiconductors by 
means of optical excitation of electrons and holes P, |3, 
which consecutively relax through mutual scatterings !3| 
or interactions with the lattice 4]. While in bulk semi- 
conductors the scattering times are completely deter- 
mined by the genuine solid-state interactions, appropri- 
ate tailoring of the carrier density of states in lower- 
dimensional semiconductor structures makes possible a 
modification of the scattering characteristics. An ex- 
treme case is found in semiconductor quantum dots |^,|^, 
often referred to as artificial atoms, where the optical 
excitation of an electron-hole pair in the state of lowest 
energy causes the deformation of the surrounding lat- 
tice but relaxation is completely inhibited because of the 
atomic-like carrier density of states. In coherent optical 
spectroscopy |i ^] , which is sensitive to the optically in- 
duced coherence, this partial transfer of quantum coher- 
ence from the electron-hole state to the lattice degrees of 
freedom, i.e., phonons, results in dephasing 0- The cou- 
pled dot-phonon system is conveniently described within 
the independent Boson model 8] , which is exactly solv- 
able for laser pulses (5-likc in time ,fl, IQsj^ whereas ap- 
proximate description schemes have to be employed for 
laser pulses of finite duration. This has recently been ac- 
complished within a density-matrix approach by Forstner 
et al. |l2l |. who reported a surprisingly large impact of 
phonon-assisted dephasing on the coherent optical re- 
sponse. Apparently, this constitutes a serious drawback 
for all-optical quantum information applications in quan- 
tum d ots, which have recently received considerable in- 
terest lilElElIilS 

It should, however, be noted that, contrary to other 
decoherence channels in solids where the system's wave- 
function acquires an uncontrollable phase through envi- 
ronment coupling in the independent Boson model 
the loss of phase coherence is due to the coupling of the 
electron-hole state to an ensemble of harmonic oscilla- 
tors which all evolve with a coherent time evolution but 
different phase. This results in destructive interference 
and dephasing, and thus spoils the direct applicability 
of coherent carrier control. On the other hand, the co- 



herent nature of the state- vector evolution suggests that 
more refined control strategies might allow to suppress 
dephasing losses. To address this problem, in this letter 
we examine phonon-assisted dephasing within the frame- 
work of optimal control 19^, 2CL, 21] aiming at a most 
efficient control strategy to channel the system's wave- 
function through a sequence of given states. We find 
that appropriate tailoring of laser pulses allows to pro- 
mote the system from the ground state through a se- 
quence of excited states back to the ground state with- 
out suffering significant dephasing losses. Despite the 
widespread use of the independent Boson model, e.g., for 
the description of optical properties of localized states in 
solids or MoBbauer spectroscopy, to our best knowledge 
no such control strategy for suppression of environment 
losses has hitherto been reported in the literature. Our 
results thus not only reinforce quantum dots as viable 
candidates for quantum- information processing devices, 
but also suggest related schemes for the control of im- 
purity or defect states in solids, e.g., in the context of 
ultraslow light propagation |22l | or creation of multipar- 
ticle entanglement '2^. 

In our theoretical approach we follow Refs. 0, 0, 0| 
and start with the usual independent Boson Hamiltonian. 
We describe the dot states in terms of a generic two-level 
system, with ground state and excited state x, assum- 
ing a negligible contribution of excited exciton states due 
to the typically large energy splittings of several tens of 
meV 01 and of biexcitons, which, in optical experiments, 
can be achieved through appropriate polarization filter- 
ing. This two-level system is coupled to a reservoir of 
harmonic oscillators such that the interaction only oc- 
curs when the system is in the upper state |^: 



H = ^gx{ax + a{)\x){x\+'^u}xa{ax 

A A 

- ^{n\x){o\ + n*\o){x\). (1) 

Here, the bosonic degrees of freedom A with energy ux are 
described by the bosonic field operators ax and a\, and 
gx is the coupling constant between x and_ 
we have neglected in accordance to Refs. 
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FIG. 1: Results of our calculations with Gaussian 2n 
(dashed lines) and optimal-control (solid lines) laser pulses 
and for: (a,d) zero temperature; (b,e) zero temperature and 
an electron-phonon coupling enhanced by a factor of \/3; and 
(c,e) finite temperature T = Wc ~ 10 K. The upper graphs 
(a-c) show |n(t)|, the lower ones (d-f) the time evolution of 
U3(t), and the insets the trajectories of the Bloch vector u{t). 
For the Gaussian 27r-pulse Rabi flopping occurs but is damped 
due to electron-phonon interactions. For the optimal control 
dephasing losses are completely suppressed, and the the sys- 
tem passes through the desired states of 63 at time zero and 
-63 at T. 



FIG. 2: Time evolution of the different phonon occupancies 
in the excited state J2\ I {^\'''>'\''Pi^-)) P S{uj-'UJx) as a function 
of phonon energy ui and for the (a) Gaussian and (b) optimal 
control shown in Fig. 1(b). 



higher-order phonon processes, e.g., anharmonic decay 
or point-defect scatterings, which is a well controlled ap- 
proximation in view of the long calculated 24] and mea- 
sured 25, 26] acoustic-phonon lifetimes. We finally de- 
scribe the light-matter coupling within the usual dipole 
and rotating-wave approximations with Q = fiS the 
Rabi frequency, the dipole moment of the two-level 
system, and S the envelope of the external laser con- 
trol. Similar to Ref. we consider a spherical dot 
model and acoustic deformation potential interactions 
9q — (9/2/32)2 (i^e — ^h) exp(— (7^(T^/4) as the only cou- 
pling mechanism, with q the phonon wavevector, p the 
mass density, c the longitudinal sound velocity, De and 
Dh the deformation potentials for electrons and holes, 
respectively, and a the carrier localization length. We 
assume a — b nm and use material parameters represen- 
tative for GaAs |0] . In the following length is measured 
in units of cr, energy in units of = c/a 0.7 meV, and 
time in units of uj^^ ^ 1 ps. 

Zero temperature. — At zero temperature the solution 
of the independent boson model reduces to that of 
the Schrodinger equation i]'0) ~ H\^), which we numeri- 
cally solve by expanding in the basis of ]0), \x), a\\0), 
and aj,la;). The control Q{t) is assumed to have a Gaus- 
sian envelope with a full- width of half maximum of 5 uj~^ 
and area 27r, and to be tuned to the polaron ground state 
frequency A = —J2xi9^)'^/'^>^ 111- To describe the two- 
level system we introduce the Bloch vector u = {a) , with 
ai the usual Pauli matrices [i^l, where ui and U2 account 
for the real and imaginary part of the interband polariza- 
tion, respectively, and U3 gives the population difference 
between and x. While in absence of phonon coupling 
n(t) would simply rotate u from —63 through a series 
of excited states back to —63, the phonon coupling of 
Eq. entangles the two-level system with the lattice 
degrees of freedom and leads to dephasing. This can be 
seen in Fig. 1(d) which shows that after the action of 



the Gaussian 27r-pulse the Bloch vector remains in an 
excited state. From the inset of Fig. 1(d), and even more 
clearly from Fig. 1(e) showing results for a phonon cou- 
pling enhanced by a factor of -s/S, we observe that the 
final deviation of u from the initial state —63 is not due 
to an incomplete rotation of u but to a loss of norm of 
]|u]j, i.e., the system has suffered dephasing losses. 

In the following we shall address the question whether 
such losses are inherent to the system under investiga- 
tion or can be suppressed by more sophisticated con- 
trol strategies. To this end, we quantify the objec- 
tive of the control through the cost function J(tt, 57) — 
lU-T^t m\u{t)-e^\'' + \u{T) + e:i\' + a j^^dt\n{t)\-'), 
with [3 a Gaussian centered at time zero with a narrow 
half-width of full maximum of {).\ijJ~^ and a = 10~^ a 
small constant, i.e, we are seeking for solutions where u 
passes through 63 at time zero and goes back to —63 
at T, and the last term in J{u,U) accounts for the lim- 
ited laser resources and is needed to make the optimal- 
control problem well posed. To determine the optimal 
control we have to minimize J(u, fi) subject to the con- 
straint that -0 fulfills Schrodinger's equation with the ini- 
tial condition]^ (—T)) = ]0), which is done in accordance 
to Refs. [23, 13 by introducing Lagrange multipliers for 
the constraints and utilizing that the Lagrange functional 
admits a stationary point at the solution. The resulting 
optimality system is solved b y u se of a gradient-related 
method including linesearch [2l|, which requires to in- 
tegrate the state equation for "ij) forwards in time and 
an adjoint equation for the dual variable ^ backwards in 
time. The solutions "0 and ip provide us with an improved 
57(i), and allow for an iterative minimization of J{u,ip) 
that converges after typically a few hundred iterations 
to the optimal control. Results of our optimal-control 
calculations are shown in Fig. 1. Most remarkably, we 
can indeed obtain a control field for which u{t) passes 
through the desired states of 63 at time zero and —63 
at T. Thus, appropriate pulse shaping allows to fully 
control the two-level system even in presence of phonon 
couplings. 

We emphasize that, contrary to the quantum "bang- 
bang" control reported by Viola and Lloyd [2I] where 
the system is constantly flipped to suppress decoherence, 
in Fig. 1 the Bloch vector is rotated slowly such that all 
oscillators can follow the system almost adiabatically ir- 
respective of Lij\. This is shown clearly in Fig. 2 where we 
plot the time evolution of the different phonon occupan- 
cies in the excited state '^^^ \ {x\a\\^l){i)) 5{uj — ijJ\). A 
comparison of the Gaussian and optimal control strate- 
gies of Figs. 2(a) and (b) reveals that only in the latter 
case all low-energy phonon modes can follow u{t) and 
become de-excited at the end of the pulse. In a sense, 
this finding is reminiscent of self-induced transparency 
I27I Efll Isfil ] where an intense laser pulse propagating in a 
medium of inhomogeneously broadened two-level systems 
acquires a pulse shape which drives all systems from the 
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groundstate through a sequence of excited states back to 
the groundstate irrespective of their detuning. 

Finite temperature. — We next address the question 
whether our findings prevail in case of finite temperature. 
Similar to Ref. |l2l |. we use a density-matrix approach 
with u = (cr) the Bloch vector, s\ = {a\) the coherent 
phonon amplitude, and u\ = {(T{a\ — s\)) the phonon- 
assisted density matrix as dynamic variables. Their time 
evolution is governed by |8 Ij : 



with the initial conditions u{—T) = —63, s\{—T) = 0, 
and ma(-T) = 0. Here, nx = {{a{ - sl){ax - sa)) is 
the occupation of incoherent phonons, which we shall 
approximate by a Bose-Einstein distribution [l^ . and 
we have introduced the abbreviation ft = — 5fe(f2) ei + 
3m(il) 62 + 2 3fJe J2x 5^ '^^ ^3- again use the method 
of Lagrange multipliers to minimize J(tt, iT) subject to 
Eq. and obtain the adjoint equations 



u = Q X u + 2 g\ 63 X jRe(MA) 

A 

SA = -it^A SA - 5 5a (1 + M3) 

iix = ft X ux — iiux ux 

+ 5a ("A + 5) 63 X It + I {U3 14-63) 



(2a) 
(2b) 

(2c) 



ii = n X u + ^gxi{nx + ^)e3 ^^e{ux) + ^"^mUsx - uux) 63 - usuxj] + f3it){u ~ e^) 



Sx = -it^ASA + '■ 



2gx({ux u) es + ne^{u*y x uy) 63 



(3a) 
(3b) 



Ma = O X Ma - icjA^A + 2 5a 63 X It , 



(3c) 



with terminal conditions m(T) = —63 — m(T), sa(T) — 0, 
and ma(T) = 0. Eqs. (0) and Q together with 

il^ - ^(m X m) + 3fie ^(m^ X ma)^ (ei - ^62) (4) 

form the optimality system. Similar to Ref. l2ll| it is 
solved iteratively through integration of Eq. ^ forwards 
and Eq. Q backwards in time, and computing an im- 
proved control by use of Eq. Q . 

At zero temperature, i.e., for nx = 0, the solutions of 
Eq. (|2Jl almost coincide with the corresponding wavefunc- 
tion results shown in Figs. l(d,e). With increasing tem- 
perature it is still possible to find optimal-control fields 
where dephasing losses are strongly suppressed, but the 
control strategy gradually changes to that reported in 
Fig. 1(f) for a temperature of T = ^ 10 K. Here, 
the transition time of u from —63 to 63 is minimized and 
consequently dephasing losses, which only occur for finite 
values of ui and M2, are suppressed. However, as com- 
pared to the low-temperature case this control strategy 
is much more trivial, and could have been guessed from 
more simple considerations. 

In conclusion, we have studied the phonon-induced 
dephasing dynamics in optically excited semiconductor 



quantum dots within the frameworks of the indepen- 
dent Boson model and optimal control. We have shown 
that appropriate tailoring of laser pulses allows to con- 
trol the dot states without suffering significant dephasing 
losses, not only at the lowest but, though exceedingly dif- 
ficult, also at elevated temperatures. The requirements 
for such laser-pulse shaping are well within the possi- 
bilities of presentday technology To highlight the 
applicability of quantum control, in this work we have 
focused on laser pulses with durations of a few picosec- 
onds where the effects of dephasing losses are most pro- 
nounced. For other control objectives it might be advan- 
tageous to use shorter or longer laser pulses, for which 
control becomes substantially simplified, or to rely on 
more advanced control strategies such as, e.g., stimu- 
lated Raman adiabatic passage 0, [s^- Besides their 
importance for future quantum-information processing 
applications, our findings might be also useful to address 
more fundamental questions regarding the nature of scat- 
terings in solids. Contrary to higher-dimensional semi- 
conductors, where scatterings occur on a sub-picosecond 
timescale 0, Qi S within the present scheme the 
atomic-like density of states leads to a significant slow- 
down of electron-phonon scatterings and allows their 
manipulation through an external control. This situa- 
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tion resembles that of controlled collisions of optically 
trapped atoms which are used to establish multi-particle 
entanglement In a similar fashion, one might envi- 
sion that controlled electron-phonon scatterings in semi- 
conductors, i.e., genuine solid-state interaction channels, 
might open a way for tunable interdot interactions. 
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sions. This work has been supported in part by the Fonds 
zur Forderung der wissenschaftlichen Forschung (FWF) 
under SRC 03 "Optimization and Control" and project 
P15752-N08. 
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